A Tutorial of the Poisson Random Field Model in Population Genetics

Population genetics is the study of allele frequency changes driven by various evolutionary forces such as mutation, natural selection, and random genetic drift. Although natural selection is widely recognized as a bona-fide phenomenon, the extent to which it drives evolution continues to remain unclear and controversial. Various qualitative techniques, or so-called “tests of neutrality”, have been introduced to detect signatures of natural selection. A decade and a half ago, Stanley Sawyer and Daniel Hartl provided a mathematical framework, referred to as the Poisson random field (PRF), with which to determine quantitatively the intensity of selection on a particular gene or genomic region. The recent availability of large-scale genetic polymorphism data has sparked widespread interest in genome-wide investigations of natural selection. To that end, the original PRF model is of particular interest for geneticists and evolutionary genomicists. In this article, we will provide a tutorial of the mathematical derivation of the original Sawyer and Hartl PRF model

Illuminating microRNA Transcription from the Epigenome

Cellular gene expression is governed by a complex, multi-faceted network of regulatory interactions. In the last decade, microRNAs (miRNAs) have emerged as critical components of this network. miRNAs are small, non-coding RNA molecules that serve as post-transcriptional regulators of gene expression. Although there has been substantive progress in our understanding of miRNA-mediated gene regulation, the mechanisms that control the expression of the miRNAs themselves are less well understood. Identifying the factors that control miRNA expression will be critical for further characterizing miRNA function in normal physiology and pathobiology. We describe recent progress in the efforts to map genomic regions that control miRNA transcription (such as promoters). In particular, we highlight the utility of large-scale “-omic” data, such as those made available by the ENCODE and the NIH Roadmap Epigenomics consortiums, for the discovery of transcriptional control elements that govern miRNA expression. Finally, we discuss how integrative analysis of complementary genetic datasets, such as the NHGRI Genome Wide Association Studies Catalog, can predict novel roles for transcriptional mis-regulation of miRNAs in complex disease etiology

Needles in the genetic haystack of lipid disorders: single nucleotide polymorphisms in the microRNA regulome

In recent years, microRNAs (miRNA) have emerged as important posttranscriptional regulators of gene expression in a wide variety of biological pathways. Since the discovery of the liver-specific miRNA-122 (miR-122) and its critical role in hepatic function, numerous additional miRNAs have been implicated in lipid metabolism. It is now apparent that lipid homeostasis is governed in part by an intricate web of miRNA activity. miRNAs are thought to confer robustness against environmental changes, such as diet modifications. Therefore, naturally occurring genetic variation that perturbs miRNA expression and/or function is likely to contribute to interindividual variability in lipid phenotypes. Although the field is still in its infancy, this review describes the growing evidence for miRNA-related genetic variation as etiological factors in lipid disorders. Specific examples, including a variant in a miRNA transcriptional control element that leads to dyslipidemia as well as a variant in a miRNA target site that modulates the effect of diet on plasma lipid levels, are discussed. Finally, the utility of recent systems genetics approaches to uncover hidden miRNA-related genetic associations with lipid disorders are considered, thereby illuminating the needles in the genetic haystack

Arsenic Exposure and Type 2 Diabetes: MicroRNAs as Mechanistic Links?

PURPOSE OF REVIEW: The goal of this review is to delineate the following: (1) the primary means of inorganic arsenic (iAs) exposure for human populations, (2) the adverse public health outcomes associated with chronic iAs exposure, (3) the pathophysiological connection between arsenic and type 2 diabetes (T2D), and (4) the incipient evidence for microRNAs as candidate mechanistic links between iAs exposure and T2D. RECENT FINDINGS: Exposure to iAs in animal models has been associated with the dysfunction of several different cell types and tissues, including liver and pancreatic islets. Many microRNAs that have been identified as responsive to iAs exposure under in vitro and/or in vivo conditions have also been shown in independent studies to regulate processes that underlie T2D etiology, such as glucose-stimulated insulin secretion from pancreatic beta cells. Defects in insulin secretion could be, in part, associated with aberrant microRNA expression and activity. Additional in vivo studies need to be performed with standardized concentrations and durations of arsenic exposure in order to evaluate rigorously microRNAs as molecular drivers of iAs-associated diabetes

Discovery of active enhancers through bidirectional expression of short transcripts

Abstract Background Long-range regulatory elements, such as enhancers, exert substantial control over tissue-specific gene expression patterns. Genome-wide discovery of functional enhancers in different cell types is important for our understanding of genome function as well as human disease etiology. Results In this study, we developed an in silico approach to model the previously reported phenomenon of transcriptional pausing, accompanied by divergent transcription, at active promoters. We then used this model for large-scale prediction of non-promoter-associated bidirectional expression of short transcripts. Our predictions were significantly enriched for DNase hypersensitive sites, histone H3 lysine 27 acetylation (H3K27ac), and other chromatin marks associated with active rather than poised or repressed enhancers. We also detected modest bidirectional expression at binding sites of the CCCTC-factor (CTCF) genome-wide, particularly those that overlap H3K27ac. Conclusions Our findings indicate that the signature of bidirectional expression of short transcripts, learned from promoter-proximal transcriptional pausing, can be used to predict active long-range regulatory elements genome-wide, likely due in part to specific association of RNA polymerase with enhancer regions

miRGen: a database for the study of animal microRNA genomic organization and function

miRGen is an integrated database of (i) positional relationships between animal miRNAs and genomic annotation sets and (ii) animal miRNA targets according to combinations of widely used target prediction programs. A major goal of the database is the study of the relationship between miRNA genomic organization and miRNA function. This is made possible by three integrated and user friendly interfaces. The Genomics interface allows the user to explore where whole-genome collections of miRNAs are located with respect to UCSC genome browser annotation sets such as Known Genes, Refseq Genes, Genscan predicted genes, CpG islands and pseudogenes. These miRNAs are connected through the Targets interface to their experimentally supported target genes from TarBase, as well as computationally predicted target genes from optimized intersections and unions of several widely used mammalian target prediction programs. Finally, the Clusters interface provides predicted miRNA clusters at any given inter-miRNA distance and provides specific functional information on the targets of miRNAs within each cluster. All of these unique features of miRGen are designed to facilitate investigations into miRNA genomic organization, co-transcription and targeting. miRGen can be freely accessed at

tDRmapper: challenges and solutions to mapping, naming, and quantifying tRNA-derived RNAs from human small RNA-sequencing data

Abstract Background Small RNA-sequencing has revealed the diversity and high abundance of small RNAs derived from tRNAs, referred to as tRNA-derived RNAs. However, at present, there is no standardized nomenclature and there are no methods for accurate annotation and quantification of these small RNAs. tRNA-derived RNAs have unique features that limit the utility of conventional alignment tools and quantification methods. Results We describe here the challenges of mapping, naming, and quantifying tRNA-derived RNAs and present a novel method that addresses them, called tDRmapper. We then use tDRmapper to perform a comparative analysis of tRNA-derived RNA profiles across different human cell types and diseases. We found that (1) tRNA-derived RNA profiles can differ dramatically across different cell types and disease states, (2) that positions and types of chemical modifications of tRNA-derived RNAs vary by cell type and disease, and (3) that entirely different tRNA-derived RNA species can be produced from the same parental tRNA depending on the cell type. Conclusion tDRmappernot only provides a standardized nomenclature and quantification scheme, but also includes graphical visualization that facilitates the discovery of novel tRNA and tRNA-derived RNA biology

Environmental contaminants and microRNA regulation: Transcription factors as regulators of toxicant-altered microRNA expression

MicroRNAs (miRNAs) regulate gene expression by binding mRNA transcripts and inhibiting translation and/or inducing degradation of the associated transcripts. Expression levels of miRNAs have been shown to be altered in response to environmental toxicants, thus impacting cellular function and influencing disease risk. Transcription factors (TFs) are known to be altered in response to environmental toxicants and play a critical role in the regulation of miRNA expression. To date, environmentally-responsive TFs that are important for regulating miRNAs remain understudied. In a state-of-the-art analysis, we utilized in silico bioinformatic analysis to characterize potential transcriptional regulators of environmentally-responsive miRNAs. Using the miRStart database, genomic sequences of promoter regions for all available human miRNAs (n=847) were identified and promoter regions were defined as −1000/+500 base pairs from the transcription start site. Subsequently, the promoter region sequences of environmentally-responsive miRNAs (n=128) were analyzed using enrichment analysis to determine overrepresented TF binding sites (TFBS). While most (56/73) TFs differed across environmental contaminants, a set of 17 TFs was enriched for promoter binding among miRNAs responsive to numerous environmental contaminants. Of these, one TF was common to miRNAs altered by the majority of environmental contaminants, namely SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin, subfamily A, member 3 (SMARCA3). These identified TFs represent candidate common transcriptional regulators of miRNAs perturbed by environmental toxicants

Prioritization of Genetic Variants in the microRNA Regulome as Functional Candidates in Genome-Wide Association Studies

Comprehensive analyses of results from genome-wide association studies (GWAS) have demonstrated that complex disease/trait-associated loci are enriched in gene regulatory regions of the genome. The search for causal regulatory variation has focused primarily on transcriptional elements, such as promoters and enhancers. microRNAs (miRNAs) are now widely appreciated as critical posttranscriptional regulators of gene expression and are thought to impart stability to biological systems. Naturally occurring genetic variation in the miRNA regulome is likely an important contributor to phenotypic variation in the human population. However, the extent to which polymorphic miRNA-mediated gene regulation underlies GWAS signals remains unclear. In this study, we have developed the most comprehensive bioinformatic analysis pipeline to date for cataloging and prioritizing variants in the miRNA regulome as functional candidates in GWAS. We highlight specific findings, including a variant in the promoter of the miRNA let-7 that may contribute to human height variation. We also provide a discussion of how our approach can be expanded in the future. Overall, we believe that the results of this study will be valuable for researchers interested in determining whether GWAS signals implicate the miRNA regulome in their disease/trait of interest

Long Noncoding RNA Moderates MicroRNA Activity to Maintain Self-Renewal in Embryonic Stem Cells

Of the thousands of long noncoding RNAs expressed in embryonic stem cells (ESCs), few have known roles and fewer have been functionally implicated in the regulation of self-renewal and pluripotency, or the reprogramming of somatic cells to the pluripotent state. In ESCs, Cyrano is a stably expressed long intergenic noncoding RNA with no previously assigned role. We demonstrate that Cyrano contributes to ESC maintenance, as its depletion results in the loss of hallmarks of self-renewal. Delineation of Cyrano's network through transcriptomics revealed widespread effects on signaling pathways and gene expression networks that contribute to ESC maintenance. Cyrano shares unique sequence complementarity with the differentiation-associated microRNA, mir-7, and mir-7 overexpression reduces expression of a key self-renewal factor to a similar extent as Cyrano knockdown. This suggests that Cyrano functions to restrain the action of mir-7. Altogether, we provide a view into the multifaceted function of Cyrano in ESC maintenance